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Using overlap as well as Wilson fermions, we have computed the one-loop renor- 
malization factors of ten non-singlet operators which measure the third moment of 
quark momentum and helicity distributions (the lowest two having been computed 
in a previous paper) , as well as the lowest three moments of the 52 structure function 
and the lowest two non-trivial moments of the h\ transversity structure function 
(plus the tensor charge). These factors are needed to extract physical observables 
from Monte Carlo simulations of the corresponding matrix elements. 



An exact chiral symmetry is maintained in our calculations with overlap fermions, 
and its most important consequence here is that the operators measuring §2 do not 
show any of the power-divergent mixings with operators of lower dimension which 
are present in the Wilson case. Many of our results for Wilson fermions are also 
new; for the remaining ones, we agree with the literature except in one case. The 
computations have been carried out using the symbolic language FORM, in a general 
covariant gauge, which turns out also to be useful in checking the gauge-invariance 
of the final results. 
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1 Introduction 



Overlap fermions [1] have emerged in recent years as one of the most promising 
formulations for simulating on the lattice theories that possess an exact chiral 
symmetry. The overlap-Dirac operator proposed by Neuberger [2] is one of the 
solutions of the Ginsparg- Wilson relation [3] 

l5 D + D l5 = a-D l5 D, (1) 

P 



and allows the realization of an exact chiral symmetry also at non-zero lattice 
spacing [4] without giving up other important symmetries. Although simu- 
lations with the Neuberger operator look computationally demanding when 
compared with Wilson fermions, progress is under way in simulating overlap 
fermions and in improving the efficiencies of their simulation algorithms [5-8]. 
Results of some recent simulations show remarkable evidence of the signatures 
of the exact chiral symmetry. For example, Ref. [7] provides strong numerical 
evidence that the pion mass approaches zero in the limit mo a — > without ad- 
ditive renormalization, in the two-dimensional Schwinger model, and progress 
is on the way for QCD; in Ref. [8] it is shown how the Gell-Mann-Oakes- 
Renner relation, which is a test of chiral symmetry at finite lattice spacing, is 
satisfied to better than 1% down to quark masses as small as m§a = 0.006. 
There is some confidence that overlap fermions will at the end turn out to 
be not too computationally demanding, and in any case one has to take into 
account all the advantages of having an exact chiral symmetry as well as the 
rapid progress in computer technology that is taking place these days. 

In this paper we present the computation of the renormalization factors of 
many operators that measure moments of various structure functions, using 
overlap fermions. The operators that we consider include all three parton dis- 
tributions that characterize the quarks in the nucleon and provide a complete 
description of quark momentum and spin at leading twist: the momentum 
distribution q(x, Q 2 ), the helicity distribution Aq(x, Q 2 ) and the transversity 
distribution Sq(x, Q 2 ). We also study the g 2 structure function, which receives 
contributions from twist-3 operators. For each one of these structure func- 
tions, we have included in our study all the moments that can be measured by 
operators for which all tensor indices are distinct. On the lattice with overlap 
fermions, all these operators are then multiplicatively renormalized. 

One of the advantages of overlap fermions is that lattice renormalization does 
not induce any mixings with operators of the wrong chirality, as instead it hap- 
pens with Wilson fermions. While this issue is of the utmost importance in 
lattice simulations concerning weak decays and quantities like the CP-violation 
parameter e'/e [9], it is also relevant here for some of the operators that mea- 



sure the #2 structure function, which for Wilson fermions mix with coefficients 
diverging as l/o, while for overlap fermions they are multiplicatively renor- 
malized (see Sect. 5). Another advantage of overlap fermions is the reduction 
of the number of independent renormalization factors in a given physical situa- 
tion, as can be seen here with the operators measuring unpolarized and polar- 
ized structure functions which differ by a 75 matrix. Significant computational 
gains are achieved in the improvement of the action (which is not needed at 
all) and of the operators, as the construction of improved operators and the 
calculation of their renormalization constants is much simpler for overlap than 
for Wilson fermions (see Sect. 4). We think that this is really an important 
point for the operators considered in this paper, as the full 0(a) improvement 
of DIS operators with Wilson fermions involves a significant amount of careful 
and cumbersome calculations. 

Among the other solutions of the Ginsparg- Wilson relation, a popular one is 
given by domain-wall fermions, where however the decoupling of the chiral 
modes is achieved only in the limit in which the number of points in the fifth 
additional dimension goes to infinity. One of the most attractive features of 
overlap fermions is that, on the contrary, chiral symmetry is fully preserved 
for any finite volume of the lattice. As the authors of Ref. [10] state, "in 
practical applications, it is our present experience that it is easier to control 
chiral symmetry violations with Neuberger's operator" . 

The renormalization with overlap fermions of the lowest two moments of the 
quark momentum and helicity distributions has been already carried out in 
Ref. [11], and we refer to that paper for the general framework and conventions, 
as well as for the Feynman rules of the overlap theory that we use here. We 
recall only that the explicit form of the overlap-Dirac operator is 
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D N = -p 1 + -== , (2) 



where 



X = D w - -p, (3) 

a 



with < p < 2r, in terms of the usual Wilson-Dirac operator 

1 



D W = - 



7m(V; + V„) - arV;V, , (4) 



1 

V M V( ;r ) = _ U(x, p)ip(x + aft) — ip(x) . (5) 



a 



Since the renormalization constants for many of the operators considered in 
this paper had not even been computed with Wilson fermions, we have cal- 
culated them for all the operators considered here with Wilson fermions too, 
and in this way we have also checked some old results in the literature, finding 
a discrepancy in one case. 

This paper is organized as follows: in Sect. 2 we introduce the various oper- 
ators of which we have computed the renormalization constants, in Sect. 3 
their renormalization on the lattice is discussed and some details about the 
perturbative calculations are given, in Sect. 4 we discuss the advantages of 
using overlap over Wilson fermions with regard to the improvement, and in 
Sect. 5 we present the complete results. In the appendices we give the results 
for the quark self-energy, both for overlap and for Wilson fermions, and for 
the individual proper diagrams. 



2 Moments of structure functions 



We have considered in this paper the lowest moments of a few structure func- 
tions which give a complete description of quark momentum and spin at lead- 
ing twist: the momentum distribution q(x,Q 2 ), measured by the F\ and Fi 
unpolarized structure functions, the helicity distribution Aq(x, Q 2 ), measured 
by the gi structure function, and the transversity distribution 5q(x, Q 2 ), mea- 
sured by the hi structure function [12,13]. We have also considered the g 2 
structure function [14] , which measures the transverse spin, is chiral even and 
contains at leading order a twist-3 piece; it is one of the most accessible higher- 
twist quantities. The hi structure function is instead chiral odd and as such 
does not arise in inclusive DIS and can be measured instead in polarized Drell- 
Yan processes. For some more detailed discussions of deep inelastic scattering 
on the lattice, see Refs. [15-20,11] and references therein. 

The moments of the various distributions are, for a given flavor Q, 



1 We have chosen our conventions in such a way that there is a manifest symmetry 
between v n and a n . The operators corresponding to v n and a n have the same number 
of covariant derivatives, corresponding chiral properties, and the same renormaliza- 
tion constant for overlap fermions [11], as well as obviously in any continuum scheme. 
The same convention was also used in Ref. [11] and can be found for example in [21], 
but in other papers like e. g. [17,19] a n has a correspondence with v n+ i. 
In all operators considered here the subscript n equals the number of symmetrized 
Lorentz indices. This means that for the operators measuring d n our convention is 
still the same as [17,19]. 



2J dxx n - 1 F 1 (x,Q 2 ) = C hn (^,g(fi))v n (fi) 
o ^ 

I Q2 

j dxx n ~ 2 F 2 (x, Q 2 ) = C 2 , n {-^, g(n)) v n (fji) 

^ 

\. \ O 2 

2 / dxx n ~ 1 g 1 (x,Q 2 ) = -E 1)n (—,g([i)^a n (ij) 

o ^ 

2 / dxx n ~ 1 g 2 (x,Q 2 ) = - — — [ E 2 , n (-^ , g{n)) d n _i(n) 



(6) 



2 n 



-£i, n (— ,£((/i))a n (/i) 



2 J dxx^hix, Q 2 ) = ^B hn (%, g(n)) tM 



where Cj jTl , E ijn and B^ n denote the appropriate Wilson coefficients in the 
OPE expansions, which can be computed in continuum perturbation theory. 
The moments of the helicity and transversity distributions are given by the 
formulae 



*n+l 



2A^q, A^g(/i) = f dxx n Aq{x, ft) 



t n+1 = 25 {n) q, 5 (n) g(/i) = / dxx n 5q{x, //), 



and the axial charge is A^w — A^°'d = gA = 1-26. 

The quark operators that correspond to the moments are 2 



(7) 



71 — 1 



o 



fJ.l---fJ.ri 



^l^D^-'-D^-tinces 



We use D=D — D, with the following lattice discretizations: 
Dp i/j(x) = — U(x, /J>)ip(x + oA) — U*{x — afi, fi)jp{x — afi) 

Zi(X - 

jp(x) Dp = — tp(x + ajl)W(x, jj) — 4>{x — ajx)U{x — aft, /j,) 



(8) 



°Mi-Mn = ( \ J ^7/ii75-D w • • • £> Mr > - traces (9) 

(. \ n-2 
2 J ^wwTs-D^, • • • IV0 - traces. 

The matrix elements of the above operators on polarized quark states are 



o H(P' «1 {w-Mn}|P, s) = 2 v n [p M1 •■•p lin - traces 

(P. S1 {M1-MIP' ^ = ~ °™ [ S (M1^2 • • -P/Xn} + traCeS ] ( 10 ) 

(P» ^[MH-mIP' «) = - rf «-i [s[ Ml P{ M2 ]P M3 • • -Pmu} + tr aces 

(P, s1°Ji{M2-Mn}IP' «) = *n-i f*[^iP{w]Pw • ' -P^} + traces 

where s M is the polarization vector of the nucleon, with s 2 = —m 2 N . 

The moments of the various distributions can be studied from first principles 
by performing lattice Monte Carlo simulations of the above matrix elements, 
which determine then the various v n , a n , d n and t n quantities. These bare 
numbers need however to be renormalized. Under renormalization, the quark 
operators corresponding to momentum and helicity distributions mix in the 
flavor singlet case with operators that measure the corresponding gluon dis- 
tributions. We consider in this paper only flavor non-singlet operators, so that 
the mixing with gluon operators is forbidden. The operators corresponding to 
the transversity distribution however do not have any gluon mixing even in 
the flavor-singlet case, as there is no gluonic transversity at leading twist, i. e. 
no chiral-odd gluon operator can be constructed. 

In a previous paper we have calculated the renormalization factors of a few 
operators which measure v 2 , i>3, a 2 and 03 [11]; here we compute the renor- 
malization constants of the operators 

O m4 ^^ 1{i D 1 D 2 D^ (11) 

V4 , e = ^{iDiDxD^ + ^ {3 D 3 D 2 D 2} tp 

-^ 1{A D 4 D 2 D 2} ^ - ^ 1{3 D 3 D l D l} ^ ) (12) 

which measure the third moment of the quark momentum distribution, and 
O a4 , d = ^ {4l5 D 1 D 2 D 3} ^ (13) 

O a4 ,e = ^7{475-D4-Dl-Dl}^ + ^ {3 -f 5 D 3 D 2 D 2} ij 

-i)-i{A^t>D A D 2 D 2} ^ - ^7{375-D3-Di-Di}^, (14) 
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which measure the third moment of the g\ quark helicity distribution. To- 
gether with the computations in Ref. [11], they complete the calculation of 
the renormalization factors of the lowest three moments of these distribu- 
tions; each continuum operator has been considered twice on the lattice by 
choosing its indices in two different ways, corresponding to two different rep- 
resentations of the hypercubic group for the same moment, which renormalize 
independently on the lattice. 

We have computed also the renormalization constants of the operators 



O dl =^l[AlbDi\i) (15) 

O d2 =i)^ [A ^D { i ] D 2 }i) (16) 

= q (^7475-01-02^ + i)lAl^>D 2 D x i) - ^Yl^iDi^ ~ i>"ii^D 2 D 4 ^ 
O d3 = 4>>y [4 >y 5 D {1] D 2 D 3} ^ (17) 

= -{i>^4l5DiD 2 D 3 ilj + ip-f4'y 5 D 1 D 3 D2^ + i)^ 4 ^D 2 D x D 3 i) 

+^7475-D 2 -D 3 Di^ + ^^ h D 3 D x D 2 i) + ^j 5 D 3 D 2 Diip 
-^^ h D 4 D 2 D 3 ^ - ^ x ^D±D 3 D 2 i) - ^ x ^D 2 D 4 D 3 i) 

-^i^D 2 D 3 D4) - i\)^ h D 3 D±D 2 i\) - ^ x ^ b D 3 D 2 D^ 

which taken together with O a2 , O a:i and O a4 determine the lowest three mo- 
ments of the g 2 structure function in Eq. (6). In the Wilson case each one of 
the O dn operators, 

ih[<T'Y5D {/ll] D lla ■ ■ ■ D^yip, (18) 

mixes, due to the breaking of chirality, with a lower-dimensional operator 
which in the continuum OPE has a mass coefficient [14], 

m q ^7[a757{Mi] jD M2 • • • £ > M „}'0, (19) 

but on the lattice this mass becomes a 1/a divergent coefficient. This mixing is 
forbidden in the overlap case, and the O dn operators are multiplicatively renor- 
malized. Thus the overlap makes a complete perturbative renormalization of 
these operators possible. 

Regarding the transversity distribution h±, we have considered the twist-two 
operators that measure the tensor charge and the lowest two non-trivial mo- 
ments, 

O tl =^4i75^ (20) 



Ot 2 =^a Hll5 D 2} ^fj (21) 

O t3 =^a 4{ll5 D 2 D 3} ^j. (22) 



The moments of the various distributions shown above are all the ones that 
can be associated with an operator which has all the indices different from 
each other. There is then no mixing due to the breaking of the Lorentz group 
to the hypercubic group. In fact, all the above operators are multiplicatively 
renormalized in the overlap case; this is true also for the particular combi- 
nations in VAte and O v ^ e in which some of the indices are equal. The Od n 
operators are not multiplicatively renormalized when using Wilson fermions, 
but their mixings are in this case due only to the breaking of chiral symmetry. 

All operators measuring higher moments of the above distributions necessarily 
have at least two equal indices. In fact, the operators which would be next 
in the ladder of moments have five Lorentz indices and therefore two of them 
have to be equal: 



V5 =^ 1{a D 1 D 1 D 2 D 3} ^j 

O a5 =^ lHl5 D 1 D l D 2 D 3} ^j (23) 

O d4 =4>'j{4'y 5 D {1] D 1 D2D 3} 4) 

U =i!(T 4{1 -f 5 D 1 D2D 3} 1p. 

This leaves open the possibility of mixing with same-dimension or (more catas- 
trophic) lower-dimension operators. We do not address here this problem, 
which is left for further studies. 



3 Perturbative renormalization 



The raw numbers extracted from Monte Carlo simulations need to be renor- 
malized to physical continuum quantities. The connection is given by 



(0 con t) 




16tt 2 



tig -*%*)]• (Of), (24) 



where 



->cont,lat\ V - ^ / r . iJfd 




(or"*) - £ I h + ~^ R T ) ■ (°H ( 25 ) 



are the continuum and lattice 1-loop expressions respectively, and the tree- 
level matrix element is the same in both cases. Their difference AR t j = R 1 ?- 1 — 
Rfj nt enters then in the gauge-invariant renormalization factors 

Zij(afx, g ) = Sij - — Q-rARijiafi), (26) 

107T Z 



which renormalize the results of Monte Carlo simulations to a continuum 
scheme. 

We choose as continuum scheme the MS scheme, since commonly the Wilson 
coefficients are computed in this scheme. On the lattice the renormalization 
condition is that the 1-loop amputated matrix elements at a certain reference 
scale fj, are equal to the corresponding bare tree-level quantities. For lattice 
matrix elements of multiplicatively renormalized operators computed between 
one-quark states, this condition reads 



(p|0 lat (//)b> 2 2 = Zo(a f ,,g )-Z; 1 (an,g )-{p\OW(a)\p)\ i00 * 

= (p\0^\a)\p)\ tr ; e 2 , (27) 

where Z^ is the wave-function renormalization, computed from the quark self- 
energy 

The 1-loop lattice matrix element of a logarithmically divergent operator O 
has the form 



(p|O( )(a)b)| 1_1OOP = (j9|O(°)(a)b)| tree x (2* 

U + ^C F ( 7o log a 2 p 2 + Vo + T + 2S 

where Vo is the finite contribution of the vertex and sails diagrams (a, b and 
c in Fig. 1), To refers to the tadpole arising from the operator (d in Fig. 1), 
and S is the finite contribution (proportional to if) of the quark self-energy of 
one leg, including the leg tadpole (e and g, or f and h, in Fig. 1). We have also 
that Cf = 2N f° r ^ e SU(N C ) gauge group. The Z Q factor for the operator 
O is then given by 

Z o (afi,g ) = 1 - ^CV (70 log a 2 /i 2 + B \ (29) 



with 

B = Vo + T + S. (30) 





(a) vertex 



(b) sail 





(c) sail 



(d) operator tadpole 





(e) leg self-energy 
(excluding tadpoles) 



(f) leg self-energy 
(excluding tadpoles) 





(g) leg tadpole 



(h) leg tadpole 



Fig. 1. The graphs contributing to the 1-loop renormalization factors of the matrix 
elements (p\0\p). The operator insertion is indicated by a circle. 
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We will call "proper" contributions the ones that do not include the self-energy 
diagrams. They correspond to the diagrams a-d in Fig. 1. 

We have used computer programs to carry out the calculations of the Feyn- 
man diagrams, given their complexities and the huge number of terms that 
are present in some cases. They utilize an ensemble of routines written in 
the symbolic manipulation language FORM, and are able to do the analytic 
calculations; Fortran codes subsequently perform the numerical integrations. 
These routines are an extension of the ones used to do calculations with the 
Wilson action in various occasions [19,22] and to compute the first two mo- 
ments of the momentum and helicity distributions with overlap fermions [11]. 
Some more details can be found in Ref. [11]. 

We would like to stress that one of the major checks of our calculations, which 
are all done in a general covariant gauge [[], is the cancellation of the gauge- 
dependent terms proportional to (1 — a) in the final numbers that connect 
the lattice to the MS scheme, as in Eq. (26). Another significant check with 
overlap fermions is that the operators V4t d and O a4 ,d have the same Z well 
within numerical integration errors, as well as V4fi and O a4<e , as expected 
from chiral symmetry. Furthermore, we verified that there is no 1/a mixing 
term in the 1-loop expressions of the operators Od x , Od 2 and Od 3 , contrary to 
what happens in the Wilson case. 

The computations of the integrals for the operators V4>e , O a4 , e and Od 3 have 
been the most demanding. We had to split the FORM outputs in several pieces 
as otherwise the Fortran programs would refuse to compile them. The running 
time for the computation of the integrals on a 60 4 grid has been in these cases 
of the order of 500 hours on a 500 MHz CPU. 



4 Improvement: overlap versus Wilson fermions 



To highlight the advantages of overlap fermions relative to Wilson fermions, 
we summarize the differences in 0(a) improvement using each method. 

Although the overlap-Dirac action possesses an exact chiral symmetry and has 
no 0(a) terms, matrix elements of operators do have order a corrections and 
therefore they need to be improved. The operators considered in this paper, 



3 The gluon propagator that we use is 



4 sin ^ sin \ 



^lE^¥l v - (1 " a) ^?¥j- (31) 
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which are all of the form O = ipOip, are made free of 0(a) corrections by 
taking [23] 

O im P = </S(l - ±-aD N ) 6 (l - ±-aD N )$. (32) 



Thus, in the overlap formulation, by rotating the spinors as in Eq. (32), op- 
erators are freed of 0(a) corrections. In Ref. [22] it has been shown that this 
recipe improves the operators to all orders of perturbation theory. Thus, full 
0(a) improvement is achieved without tuning any coefficients. Furthermore, 
the renormalization constants for the improved and unimproved operators are 
the same, so there is no need to do additional calculations, which are generally 
cumbersome, to compute the renormalization factors in the improved theory. 
What happens is that in 1-loop amplitudes a factor Dn can combine with a 
quark propagator, but since it has an a in front and (contrary to the Wilson 
case) there is no 1/a piece in the propagator, as additive mass renormalization 
is forbidden by chiral symmetry, the contribution of D^ to the renormalization 
factors is zero [24]. Thus, improved Monte Carlo simulations are performed 
with the operator (32), nevertheless the renormalization of this operator is 
equal to that of the corresponding unimproved operator. 

For Wilson fermions, instead, the improvement of operators looks more compli- 
cated and troublesome. First of all, one has to find out a complete basis which 
includes all operators which are one-dimension higher than the original one 
and have the same symmetries. Then, to get the improved renormalization 
factors, which are different from the unimproved ones, one has to compute 
the 1-loop matrix elements of each one of those operators. Lastly, one has 
to determine somehow the values of the coefficients in front of the operator 
counterterms, and this appears to be a highly non-trivial task. 

Let us consider, for example, the first moment of the quark momentum distri- 
bution, V2 = ip'-f^D^ip. The improvement counterterms are two dimension- 
five operators [22], and one basis for the improvement is given by 

1 1 

ip^Dvyip - -aici(^o) ^i>o x{ ^\D v} ,D^i> - -ac 2 (gl) ^{D^ D„}^.(33) 



The operator is fully 0(a) improved only for some particular values of the 
coefficients d(g%) = 1 + gfc? + 0(g$) and c 2 (gl) = 1 + g^ + 0(g$). How- 
ever, presently one knows only a relation between the two coefficients [22], and 
thus one of them remains unknown. Although one could determine all coef- 
ficients using some Ward Identities or physical conditions, this involves more 
effort than with overlap fermions. In addition, the coefficients are different 
for different orders in perturbation theory, and the basis itself is different for 
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each operator to be improved. It is reasonable to expect that operators which 
contain more covariant derivatives will have larger improvement bases, and 
thus many more coefficients to be determined, besides computing the 1-loop 
matrix elements for each of the operator counterterms. Presumably already 
implementing the improvement for the second moment operators with Wilson 
fermions will turn out to be a daunting task. This means that in most cases it 
will be difficult to go beyond tree-level improvement for the structure function 
operators in the Wilson case. 

In addition to improving operators, one should also note that Wilson fermions 
need corrections to their action (the Sheikholeslami-Wohlert term), while the 
overlap action is already 0(a) improved. However, this is less of a complication 
than the improvement of the operators. 

Thus, for all these reasons, from the perspective of renormalization, mixing, 
and improvement, overlap fermions offer compelling advantages relative to 
Wilson fermions. 



5 Results 



We give in this section the results for the renormalization factors of every 
operator considered in this paper, both for the overlap and the Wilson action, 
for r = 1. 

In the overlap case, we have explicitly verified that the renormalization con- 
stants of O a ^d and a4 ,e are equal to the ones of V4t d and V4>e respectively, 
which are given below. This is a consequence of the exact chiral symmetry 
which the overlap theory possesses. Another important gain coming from chi- 
ral symmetry is that the operators O^ , Od 2 and Od 3 are multiplicatively renor- 
malized, while for Wilson fermions the breaking of chiral symmetry causes a 
mixing with operators of lower dimensions, with the mixing coefficients going 
to infinity in the continuum limit. 



5. 1 Overlap fermions 



We first consider the 1-loop contributions of the proper diagrams (a-d in 
Fig. 1), which are 



^proper 5*0 /m 



i 

16tt 2 



(— + (!-«)) log aV 
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V "J 1 

V4, a 


V a=1 

V4,e 


y dl a=1 


y d2 a=1 


V. a=1 

d 3 


y t a=l 


Vt *=i 


y t a=l 


0.2 


-9.61600 


-7.98774 


52.66717 


45.13633 


41.92810 


5.38049 


-4.02560 


-7.79945 


0.3 


-9.46984 


-7.84341 


33.37857 


26.64726 


23.65026 


5.19029 


-4.00839 


-7.70244 


0.4 


-9.34521 


-7.72067 


24.03566 


17.76727 


14.91522 


5.03112 


-3.99077 


-7.61751 


0.5 


-9.23484 


-7.61223 


18.60963 


12.65715 


9.91543 


4.89210 


-3.97274 


-7.54054 


0.6 


-9.13474 


-7.51412 


15.11256 


9.39617 


6.74341 


4.76733 


-3.95427 


-7.46925 


0.7 


-9.04244 


-7.42384 


12.70136 


7.17180 


4.59344 


4.65325 


-3.93535 


-7.40221 


0.8 


-8.95626 


-7.33973 


10.95882 


5.58297 


3.06844 


4.54744 


-3.91595 


-7.33847 


0.9 


-8.87501 


-7.26060 


9.65553 


4.40974 


1.95101 


4.44823 


-3.89606 


-7.27733 


1.0 


-8.79780 


-7.18557 


8.65526 


3.52185 


1.11258 


4.35437 


-3.87564 


-7.21828 


1.1 


-8.72396 


-7.11394 


7.87218 


2.83746 


0.47256 


4.26491 


-3.85467 


-7.16091 


1.2 


-8.65295 


-7.04517 


7.24967 


2.30276 


-0.02200 


4.17912 


-3.83313 


-7.10491 


1.3 


-8.58433 


-6.97884 


6.74895 


1.88097 


-0.40718 


4.09639 


-3.81099 


-7.05001 


1.4 


-8.51774 


-6.91457 


6.34259 


1.54618 


-0.70839 


4.01623 


-3.78821 


-6.99600 


1.5 


-8.45289 


-6.85207 


6.01069 


1.27965 


-0.94394 


3.93824 


-3.76477 


-6.94269 


1.6 


-8.38952 


-6.79107 


5.73851 


1.06752 


-1.12735 


3.86206 


-3.74064 


-6.88993 


1.7 


-8.32742 


-6.73135 


5.51486 


0.89934 


-1.26882 


3.78741 


-3.71578 


-6.83758 


1.8 


-8.26640 


-6.67272 


5.33115 


0.76710 


-1.37612 


3.71403 


-3.69017 


-6.78553 



Table 1 

The Feynman-gauge constants V G a=1 for the momentum, helicity and transversity 
operators considered in this work, in the overlap theory. Note that in this case 

V a=1 



V a j l 

a4,d 



V^I 1 and V a ?= l 



V4,e 



,0 proper 
V4,e 
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U dj 
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16tt 2 



16tt 2 
16tt 2 
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■c f 



■c f 



+V V4 a f - (1 - a) 7.553824 + T m4 



^jtree 
^V4,di 



(— + (!-«)) log aV 



2„2 



+ K 



a=l 



'1- a) 8.024764 + T V 



A^tree 



(34) 



a log oV + V di a=1 - (1 - a) 7.850272 + T dl 



Otrce 
di ' 



+ (1 - a)) logo p 



,2 2 
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operator 


operator tadpole 


l-A>4 ,di *J(i4,d 


T VA , d = T a4td = 16vr 2 ( - h + (1 - a)( - ±Z - \z x + 1)) 




^VAfil U(l4,e 


T VA , e = T a4 , e = i6vr 2 ( - ^Z - -Zi + - + (1 - a)(gZ + Z x - 


i)) 


o dl 


1 1 

T dl =l67r 2 (--Zo + (l-a)-Z ) 




o d2 


T d2 = 167r 2 (-Z + (l-a)iz ) 




o d3 


T, 3 = 16 7 r 2 (-^ + (l-a)(-lz -Iz 1 + l)) 




o tl 


T tl =0 




o t2 


T t2 = 167r 2 (-iz + (l-a)iz ) 




Ot 3 


T t3 = 167T 2 (-Zo + (l-a)iz ) 





Table 2 

The operator tadpoles for the various operators, where Zq = 0.154933390231 

and Zi = 0.107781313540 •••. 



+V,,? - (1 - a) 8.369693 + T c 
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U d 3 



16tt 2 



■c f 



'd 2 



(— + (!-«)) log aV 
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d-2 
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U d 2 



* 



.otree 
U d 3 



^proper 



167T 2 



C f 



2J2 



1 — a) log a p 



+^ a=1 - (1 - a) 3.792010 + T h 



oil 
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n proper _ 
U t 2 - 


9l r 

t:^F 

16vr 2 


(2 + (l-a))logoV 






+Vt 2 a=1 - (1 - «) 6.350272 + T t2 


.otree 
ty *2 


r\ proper 
U t 3 - 


9o r 
±bir z 


(— + (l-a))logaV 






+ v t 3 a=1 ~ (1 - «) 7.036360 + T <3 


^tree 



The Feynman gauge results V a=1 for the contributions of the sails and ver- 
tices, Vo, are tabulated in Table 1 for several values of the parameter p. They 
are also given separately for each diagram in Appendix B (for p = 1). The 
remaining parts proportional to (1 — a) are instead independent of p, and 
as their analytic expressions are very complicated functions of p containing 
thousands of terms, the numerical cancellation of this dependence is a highly 
non-trivial check of our computations. Furthermore, the parts proportional to 
(1 — a) have also the same value as with the Wilson action [11]. In fact they 
depend only on the gluonic action chosen. 

The result for the tensor charge O tl = ipa^^ip turns out to be equal to the 
result for the standard tensor current ipa^ip, which has been calculated in 
Ref. [24] in Feynman gauge and in Ref. [9] in a general covariant gauge. For 
Wilson fermions, the renormalization of Of 1 can be found in [22] . 

The results for the remaining proper diagram, the operator tadpole To (di- 
agram d in Fig. 1), are shown in Table 2. To compute the renormalization 
factors, one finally adds to the proper diagrams the 1-loop amplitudes of the 
self-energy and tadpole of one leg which are proportional to ifi. Their values 
are given in Appendix A. Putting everything together, we get the expressions 
of the renormalized operators on the lattice for overlap fermions, which for 
p = 1 are: 



O 



V4 , 



o 



V4,e ' 

o dl 

o d2 
o d3 



1 - 7fW^7 lo S a V - 83.12758+ ^(1 -a) 
lb7T z \ 30 6 



157 



11 



1 - S^ C f ( ^77 log oV - 85.33588 + — (1 - n ) 
167^ \ 30 6 



1 + -|^C F • 41.20842 OX c 

1 - tJ^CV (X log oV - 58.57488 + 1(1 - a) 



9l 



1 35 5 

1 j^Cf - log aV- 73.21719 + -(1 - a) 



Otree 
v 4 ,d 



jotrcc 
i>4,e 



A-itrcc 
U d 2 



^trcc 
U d 3 



(35) 



16 



o t2 

Ot 3 



1 - 

1 

1 



167T 



C F log a V - 33.27626 + (1 - a) 



/-jtree 



-|%C F ( 3 log aV - 53.73932 + ^(1 - a) 

167T Z \ 2 



.otrce 
^2 



13 



11 



j^C F f y log a 2 /i 2 - 69.31500 + y (1 - a) 



/^)trcc 



We see that in the overlap case the operators O^, Od 2 and Od 3 are multi- 
plicatively renormalized (contrary to what happens with Wilson fermions, see 
below) . 

Another check of our calculations is that the transversity operators, which 
apart from the trivial tensor charge have never been computed before on the 
lattice, agree with the 1-loop anomalous dimension formula [13] 



lu 



3=1 J 



(36) 



which also implies that the anomalous dimensions are positive and for any 
given moment greater than the ones of the corresponding moments of i*\, F 2 
and gi. 

To complete the connection with the continuum physics as in Eq. (24), we 
need also the 1-loop results for the same matrix elements in the continuum 
MS scheme, which are given by 



/)MS 



1- 



o 



MS 



16tt 2 

A. 

16vr 2 



■c f 



■c f 
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log a fi 



2216 11 



225 
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i 2 2 2216 11., 

log a a 1 (1 
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1- 
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1- 
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/=jMS 
U t 3 - 


1- 
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.otree 
^124 
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(37) 
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The connection of overlap lattice fermions with the continuum MS is then 
given by the gauge-invariant factors 
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u v 4 ,d 
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V4,e 



,olat (overlap) 
a4,e 



log oV - 56.63044 
— log oV- 69.80979 



O 



lat (overlap) 



lat (overlap) 



logaV- 32.27626 



) 



O 



o 



lat (overlap) 
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O 



lat (overlap) 
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For /3 = 6.0, // = l/o and N c = 3 one has 
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«2 
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= 0.99559 ojf (Wilson) 


lat 
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overlap) 


= 1.10021 ojf (Wilson) 



(39) 



where we have also shown the corresponding Wilson results (see Eq. (42)). 
We remind that, although when p = 1/a the logarithms disappear from the 
renormalization factors, for a general p they are still present but they cancel 
then against the corresponding logarithms in the Wilson coefficients, as the 
moments Eq. (6) have to be independent of the renormalization scale (and all 
the scale dependency is in the logarithms). 

The renormalization constants for overlap fermions reported here look in gen- 
eral large, especially when compared to the corresponding Wilson results, as 
shown in Eq. (39). However, if one looks at the overlap results for the proper 
diagrams (see for example Appendix B), one can notice that they are not so 
much different from Wilson fermions. In general, the biggest contribution to 
the renormalization constants comes from the operator tadpoles, but as they 
are exactly the same for overlap and Wilson fermions, it is not in these dia- 
grams that the difference can be found. One instead has to look at the quark 
self-energy (see Appendix A). In the Feynman gauge, the leg self-energy in the 
overlap (for p = 1) is -37.63063, while for Wilson fermions is +11.85240; their 
difference is -49.48303, and it is remarkable how close this number comes 
to the differences between the complete finite contributions Bq (defined in 
Eqs. (29) and (30)) for the two kinds of fermions, Eq. (38) vs. Eq. (42). If one 
would consider the overlap for p = 1.9, the difference between the self-energies 
would go from -49.48303 down to -26.80898; however, the quark propagator 
becomes singular for p = 2, so simulations would likely be more expensive 
when approaching this value of p. 



5.2 Wilson fermions 



The 1-loop contributions of the proper diagrams are in the Wilson case 
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where 07^ = §[7^,7,,]. 



The operator tadpoles are the same as in the overlap case (see Table 2). Adding 
the quark self-energy (appendix A), the complete renormalization factors are 
then 
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We can see that for Wilson fermions the operators O^, Od 2 and Od 3 are not 
multiplicatively renormalized. Each of them mixes with an operator which 
is one dimension lower [17], and the corresponding mixing coefficients, which 
seem to be gauge- invariant, are proportional to r/a. Thus, these mixings would 
be zero for naive fermions, and they are akin to the S term in the Wil- 
son quark self-energy, which is responsible for the additive renormalization of 
quark masses when they are not anymore protected by chirality. In fact, the 
lower-dimensional operators above are all mass terms of the form (19) coming 
from the OPE expansions in DIS. All this hints to a connection with chirality 
in the mixings of the Od n operators too, and we have indeed shown in Eq. (35) 
that in the calculations done with overlap fermions, where chirality is con- 
served, there is no trace of this kind of terms and these operators are then 
multiplicatively renormalized. 

We can also notice from the Wilson results above that another consequence of 
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the breaking of chiral invariance is that the renormalization constants of V4} d 
and O ai ^ are not equal, and the same happens for V4fi and O a4j£ . 

The perturbative calculations with Wilson fermions of the renormalization 
constants for the operators V4t d, O a4j d, V4<e , Od 2 and O tl have been already 
done in the past in Feynman gauge. We agree with the results for V4t d and 
O a4t d in Ref. [25], but we find discrepancies with the renormalization of 0„ 4 , e 
in Ref. [17] 0. We agree also with the results for Od 2 in Ref. [17] and for O tl 
in Ref. [22]. 

The connection of Wilson fermions with the continuum MS is given by 
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(42) 



and the numbers for (3 = 6.0, fi = 1/a and N c = 3 are given in Eq. (39). We 
have not included here Od 11 Od 2 and Od 3 , as their renormalization involves 
mixing coefficients that diverge in the limit of zero lattice spacing, and as 
such it is of no use to compute the connection to MS for these operators in 
perturbation theory with Wilson fermions. 



4 In this case we do not even agree on the results for the individual proper diagrams, 
which are given in Table 5 of Ref. [17]. In that Table, the operator tadpole is given 
the value that would be appropriate for V4j d, but this cannot be true for O v 



'U4,ej 



since there are additional Wick contractions between the A^s originating from the 
equal indices. 
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6 Conclusions 



Overlap fermions are one of the most promising formulations for putting chi- 
ral fermions on the lattice and for studying long-standing problems linked 
with chirality. In this paper we have computed the renormalization constants 
of the lowest moments of various structure functions which give a complete 
description of the quark momentum and spin at leading twist. We have com- 
puted these constants also with Wilson fermions, since the renormalization of 
many of these operators had never before been computed on the lattice. Chi- 
ral symmetry plays an important role in the structure of the strong radiative 
corrections. In particular, the overlap is the only case in which it has yet been 
demonstrated that all the operators we have considered are multiplicatively 
renormalized. 

The numbers presented here are also valid in the unquenched case if one 
deals with flavor non-singlet quark operators, which do not mix with gluon 
operators and for which internal quark loops never have the chance to come to 
play at 1-loop level. However, the numbers for the transversity operators can 
be considered unquenched also for flavor singlet quark operators, since there 
are no gluon operators with the same quantum numbers. 

We remind again that all the renormalization constants presented in this pa- 
per are already fully 0(a) improved in the overlap case. In the Wilson case 
removing all order a effects would instead involve a great amount of additional 
calculations. 

With overlap fermions 1-loop corrections are substantially larger than for Wil- 
son fermions, with the primary contribution arising from quark self-energies. 
This is an important physical effect that should be understood, and may ulti- 
mately suggest an appropriate form of resummation or tadpole improvement. 
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p 


non-tadpole self-energy 


leg tadpole 


total self-energy 


0.2 


-27.511695 +11.911596 f 


-213.087934 -7.119586 £ 


-240.599629 


+4.792010 £ 


0.3 


-23.687573 +11.098129 £ 


-131.723110 -6.306119 £ 


-155.410693 


+4.792010 £ 


0.4 


-21.172454 +10.520210 £ 


-91.817537 -5.728200 £ 


-112.989991 


+4.792010 £ 


0.5 


-19.337313 +10.071356 £ 


-68.315503 -5.279346 £ 


-87.652816 


+4.792010 £ 


0.6 


-17.912921 +9.704142 £ 


-52.931363 -4.912132 £ 


-70.844284 


+4.792010 £ 


0.7 


-16.760616 +9.393275 £ 


-42.140608 -4.601265 £ 


-58.901224 


+4.792010 £ 


0.8 


-15.800204 +9.123666 £ 


-34.193597 -4.331656 £ 


-49.993801 


+4.792010 £ 


0.9 


-14.981431 +8.885590 £ 


-28.125054 -4.093580 £ 


-43.106485 


+4.792010 £ 


1.0 


-14.270881 +8.672419 £ 


-23.359746 -3.880409 £ 


-37.630627 


+4.792010 £ 


1.1 


-13.645294 +8.479438 £ 


-19.534056 -3.687428 £ 


-33.179350 


+4.792010 £ 


1.2 


-13.087876 +8.303183 £ 


-16.407174 -3.511173 £ 


-29.495050 


+4.792010 £ 


1.3 


-12.586126 +8.141044 £ 


-13.813486 -3.349034 £ 


-26.399612 


+4.792010 £ 


1.4 


-12.130497 +7.991018 £ 


-11.635482 -3.199008 £ 


-23.765979 


+4.792010 £ 


1.5 


-11.713524 +7.851554 £ 


-9.787582 -3.059544 £ 


-21.501106 


+4.792010 £ 


1.6 


-11.329238 +7.721442 £ 


-8.206069 -2.929432 £ 


-19.535307 


+4.792010 £ 


1.7 


-10.972744 +7.599750 £ 


-6.842630 -2.807740 £ 


-17.815374 


+4.792010 £ 


1.8 


-10.639905 +7.485778 £ 


-5.660084 -2.693768 £ 


-16.299989 


+4.792010 £ 


1.9 


-10.327042 +7.379023 £ 


-4.629539 -2.587013 £ 


-14.956581 


+4.792010 £ 



Table A.l 

Results for the quark self-energy with overlap fermions, where we have used the 
abbreviation £ = 1 — <x The first column (non-tadpole self-energy) refers to diagram 
e (or f) in Fig. 1, and the second column (leg tadpole) refers to diagram g (or h). 



A Quark self-energy 



In this appendix we report the results of 1-loop calculations regarding the 
contribution proportional to ift of the quark self-energy, also known as Ei. 
This corresponds to the diagrams e (or f) and g (or h) in Fig. 1, and is 
necessary for the calculation of the renormalization constants of the operators 
considered in this work. 
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For overlap fermions we have that the total self-energy is given by 



S overlap i/o /~i 

l07T 2 



a log aY + S a=l + (1 - a) 4.792010 



(A.l) 



where the Feynman-gauge finite results S a=1 are given in the last column of 
Table A.l (in which we have used the abbreviation £ = 1 — a), which also 
contains the results of the individual diagrams. For the total self-energy (but 
not for the individual diagrams) the finite part proportional to (1 — a) is 
independent of the p parameter, and even of the fermion action used. In fact 
this number comes from integrals which have only gluonic propagators, and 
is given by F - lE + 1 = 4.792009568973 • • • [9]. 

In the Wilson case, the value of the leg self-energy (including the tadpole) is 



yi Wilson jO r~t 

1 "" 16tt 2 F 



a log ay + 11.852404 + (1 - a) 4.792010 



(A.2) 



The result for the individual diagrams is the following: for diagram e (or f) is 



16tt 2 



c t 



a log a 2 p 2 - 0.380646 + (1 - a) 7.850272 



(A.3) 



while the leg tadpole, i.e. diagram g (or h), is 



qICf 



—Z — (1 — a)—Zn 
2 o v v g u 



16tt 2 



■c f 



12.233050 -(1- a) 3.058262 



,(A.4) 



where Z n = 0.154933390231 • • •. 



B Results of the individual diagrams 



In this appendix we give the results for the finite parts of the proper diagrams: 
the vertex (diagram a in Fig. 1), the sails (diagrams b plus c) and the operator 
tadpole (diagram d). We have used the abbreviation £ = I — a. Table B.l 
refers to overlap fermions (for p = 1) and Table B.2 to Wilson fermions. For 
the latter, Table B.3 also reports the numbers giving the mixing with the lower 
dimensional operators. 

Note that the sum of the (1 — a) parts of vertex and sails is independent 
of the fermion action used, and in general this is also true for the individual 
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vertex (diagram a) 


sails (diagrams b + c) 


op. tadpole (diagram d) 


V i y d 


0.70281 -0.53019 £ 


-9.50061 -7.02363 £ 


-36.69915 +4.59515 £ 


O a A,d 


0.70281 -0.53019 £ 


-9.50061 -7.02363 £ 


-36.69915 +4.59515 £ 


O v A,e 


0.61867 -0.05926 £ 


-7.80424 -7.96551 £ 


-40.51969 +5.06609 £ 


O a 4,e 


0.61867 -0.05926 £ 


-7.80424 -7.96551 £ 


-40.51969 +5.06609 £ 


o dl 


2.81906 -0.65840 £ 


5.83620 -7.19187 £ 


-12.23305 +3.05826 £ 


o d2 


1.15709 -0.49753 £ 


2.36476 -7.87216 £ 


-24.46610 +4.07768 £ 


o d3 


0.64918 -0.45310 £ 


0.46340 -8.10072 £ 


-36.69915 +4.59515 £ 


On 


4.35437 -3.79201 £ 








On 


1.46711 -1.23375 £ 


-5.34275 -5.11652 £ 


-12.23305 +3.05826 £ 


On 


0.68662 -0.54766 £ 


-7.90490 -6.48870 £ 


-24.46610 +4.07768 £ 



Table B.l 

Results of the proper diagrams for overlap fermions, where we have used the abbre- 
viation £ = 1 — a. 





vertex (diagram a) 


sails (diagrams b + c) 


op. tadpole (diagram d) 


O v A,d 


0.842048 -0.530195 £ 


-9.201715 -7.023629 £ 


-36.699150 +4.595148 £ 


O a A,d 


0.465704 -0.530195 £ 


-9.201715 -7.023629 £ 


-36.699150 +4.595148 £ 


O v A,e 


0.842455 -0.059255 £ 


-7.527094 -7.965509 £ 


-40.519687 +5.066088 £ 


O a A,e 


0.355884 -0.059255 £ 


-7.527094 -7.965509 £ 


-40.519687 +5.066088 £ 


o dl 


2.590817 -0.812814 £ 


-1.845174 -7.037458 £ 


-12.233050 +3.058263 £ 


o d2 


0.988791 -0.585340 £ 


-4.052542 -7.784353 £ 


-24.466100 +4.077683 £ 


o d3 


0.536160 -0.521292 £ 


-5.438520 -8.032532 £ 


-36.699150 +4.595148 £ 


Oa 


4.165675 -3.792010 £ 








O t 2 


0.980373 -1.233747 £ 


-5.077267 -5.116525 £ 


-12.233050 +3.058263 £ 


Ota 


0.465618 -0.547660 £ 


-7.609564 -6.488700 £ 


-24.466100 +4.077683 £ 



Table B.2 

Results of the proper diagrams for Wilson fermions, where we have used the abbre- 
viation £ = 1 — a. 
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vertex (diagram a) 


sails (diagrams b + c) 


op. tadpole (diagram d) 


o dl 


1.508200 


14.735561 





o d2 


0.180285 


4.085395 





Odx 


0.073704 


2.808116 






Table B.3 

Results for the mixing with the lower-dimensional operators for Wilson fermions. 

diagrams. An exception to this are the Od„ operators and the quark self-energy, 
and these are also the same cases for which with Wilson fermions there is a 
1-loop mixing with a l/o divergent coefficient, while with overlap fermions 
such mixings are forbidden by chiral symmetry. 
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